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Using density functional methods we calculated structural and electronic properties of bulk chloroform and 
bromoform intercalated C^q, C^q ■ 2 CHX3 (X=Cl,Br). Both compounds are narrow band insulator materials 
with a gap between valence and conduction bands larger than 1 eV. The calculated widths of the valence and 
conduction bands are 0.4-0.6 eV and 0.3-0.4 eV, respectively. The orbitals of the haloform molecules overlap 
with the n orbitals of the fullerene molecules and the p-type orbitals of halogen atoms significantly contribute 
to the valence and conduction bands of C50 ■ 2 CHX3. Charging with electrons and holes turns the systems 
to metals. Contrary to expectation, 10 to 20 % of the charge is on the haloform molecules and is thus not 
completely localized on the fullerene molecules. Calculations on different crystal structures of C60 ■ 2 CHCI3 
and C60 • 2 CHBr3 revealed that the density of states at the Fermi energy are sensitive to the orientation of the 
haloform and Cgo molecules. At a charging of three holes, which corresponds to the superconducting phase of 
pure C50 and C^q ■ 2 CHX3, the calculated density of states (DOS) at the Fermi energy increases in the sequence 
DOS(C 50 ) < DOS(C 6( , • 2 CHCI3) < DOS(C 60 ■ 2 CHBr 3 ). 
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I. INTRODUCTION 

This paper is concerned with structural and electronic prop- 
erties of solid C60 intercalated with haloform molecules, 
CHX3 (X=Cl,Br). These compounds have been studied re- 
cently by X-ray diffraction [|J ^, |^] and solid-state nuclear 
magnetic resonance (NMR) spectroscopy [Q, 01. 

These kind of fullerene based materials have attracted a lot 
of attention because of their unusual conducting and super- 
conducting properties. Solid Cgo is insulating but can be made 
conducting and even superconducting upon intercalation of al- 
kali atoms (e.g. K, Rb) between the fullerene molecules (see 
the reviews [|| |(| ^]). The alkali atoms transfer their valence 
electrons into the conduction bands of the C60 subsystem and 
the materials become metallic. 

A noticeable enhancement of the superconducting critical 
temperature, T c , of C60 was expected when electrons are re- 
moved (hole doping) because of the larger density of states 
of the valence bands. Doping of ~ 3 holes per C60 molecule 
and corresponding metalization was achieved with the use of 
field-effect doping technique and has resulted in T c = 52 K 
Such doping has also been realized recently [[)], |j~0| ] . The 
intercalation of Cgo with haloform molecules has led to an 
further increase of the critical temperature [p~T|] . Such an en- 
hancement of T c upon intercalation of molecules into a super- 
conductor has been already observed earlier jl2| ] and was dis- 
cussed theoretically Jl3[]. In this work we focus on the normal 
properties of haloform intercalated Cgo- 

X-ray powder diffraction measurements at various tem- 
peratures [jj], ||] have shown that above ~ 200 K the crys- 
tal structure of haloform intercalated C(,o is hexagonal and 
transforms to a triclinic structure below T ~ 150 K. The sto- 
ichiometric composition of these materials is C60 ■ 2 CHX3 
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(X=Cl,Br). Solid state NMR measurements on C 60 • 2 CHBr 3 
have shown that above 193 K the Cgo molecules rotate isotrop- 
ically [Q], similar to that found in pure face-centered cubic 
(fee) C60- Whereas below this temperature a rotation barrier 
of 6 kJmoU 1 was evaluated . Also the bromoform molecules 
are found to be motionally active down to 2 1 8 K [Q| . 

Dinnebier et al. [^J calculated the electronic structure 
of these materials using a two-dimensional periodic tight- 
binding formalism where only one specific lattice plane was 
considered. They omitted the haloform molecules completely 
and took only the valence electrons of the five-fold degen- 
erate h u molecular orbital (MO) of each C60 explicitely into 
account. The main conclusion of their calculations is that 
the electronic density of states (DOS) at the Fermi energy 
for ~ 3 holes per C60 molecule satisfies following sequence: 
DOS(C 60 • 2 CHBr 3 ) « DOS(C 6 o • 2 CHCI3) < DOS(C 60 ). 
Hence, the observed enhancement of T c cannot simply be ex- 
plained by an increase of the DOS at the Fermi energy when 
C60 is intercalated with haloform molecules and CHCI3 is re- 
placed by CHBr3 (vide supra). Given these observations and 
results, it is of interest to study the electronic properties of 
these materials in more detail. 

We have performed all-electron density functional (DF) cal- 
culations on several crystal structures of bulk C60 ■ 2 CHX3 
(X=Cl,Br) to predict their electronic structures and to examine 
the effect of introduced charge carriers (doping with electrons 
as well as with holes). Employed are methods that rely on 
the generalized gradient approximation (GGA) and the spin- 
unrestricted Kohn-Sham approach. Different charge doping 
levels are simulated simply by removing or adding electrons 
from the systems. Unlike previous work ||3]j our DF calcu- 
lations include all electrons of the system into the quantum 
treatment and the haloform molecules are taken into account 
explicitly. Particular attention is paid to the interaction of the 
haloform molecules with the electronic structure of the C60 
subsystem. DF calculations have also been performed on an 
fee and on a hypothetical hexagonal lattice of C60 molecules. 
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TABLE I: Crystal structures of pure Ceo an d C^o ■ 2 CHX 3 
(X=Cl,Br) reported in the literature and used in this work (unit cell 
lengths in A, angles in degrees). Note, for Cgo the lattice parameters 
of the primitive cell is given. 





pure C 60 


Qo ■ 2 


CHX 3 (X=Cl,Br) 




temp, range [K] 


170 


>200 




< 150 




Bravais lattice 


fee 


hexagonal 


triclinic 




Space groups 


Fm3 


P3, Pi, 


PI 


Pi 




acronym 


iTT3i 


HT a ) 




LT b ) 




Ref. 


i 




i 




lattice parameters 
















X=C1 


X=Br 


X=C1 


X=Br 


a 


9.949 


10.08 


10.212 


9.8361 


9.8982 


b 


9.949 


10.08 


10.212 


10.0906 


10.3386 


c 


9.949 


10.11 


10.209 


9.8179 


9.8993 


a 


60.0 


60.0 


60.0 


101.363 


100.951 


P 


60.0 


60.0 


60.0 


116.457 


115.920 


Y 


60.0 


120.0 


120.0 


79.783 


78.202 



High temperature structure. b ) Low temperature structure. 



The latter system is formed by artificially removing the halo- 
form molecules from C60 ■ 2 CHX3. These calculations allow 
to examine changes in the electronic structure of the C60 sub- 
system (i) upon intercalation and (ii) due to electronic con- 
tributions of the haloform molecules separately. Note, that 
all structures and electronic properties are obtained in the ab- 
sence of an applied electric field. 

The structure of the article is as follows: Sec. |j] summa- 
rizes the experimental determination of the crystal structure 
of C 6 o • 2 CHX 3 (X=Cl,Br) Sec. pl| gives details about the 
computational method. In Sec. 1V A and Sec. IV B we discuss 
the calculated crystal and e lectro nic structures of the un-doped 
systems, respectively. Sec . IV C then discusses the changes of 
the electronic structure due to hole- and electron-doping of 
C60 • 2 CHX3. Finally, a concluding discussion is found in 
Sec.0 



II. CRYSTAL STRUCTURES FROM X-RAY 
DIFFRACTION MEASUREMENTS 

The first crystal structures of CHCI3- and CHBr3- interca- 
lated Cgo reported in the literature were published by Jansen 
and Waidmann [[!]]. An X-ray powder diffraction analysis 
at room temperature yielded for both compounds a primi- 
tive hexagonal lattice and an average symmetry of P6/mmm. 
However, due to orientational disorder of the molecules the 
atomic coordinates could not be resolved. Throughout this 
paper this structure is labeled with HT (high temperature) and 
the observed lattice parameters are summarized in Table ffl. 

Collins et al. examined the structure of Cgo ■ 2 CHBr3 at 
room temperature using X-ray powder diffraction as well and 
found a hexagonal lattice with unit cell parameters consistent 
with those for the HT structure of Table |. They calculated 
the atomic positions using molecular modeling techniques by 
placing a C60 molecule at the origin of the hexagonal unit cell 



and inserting the CHBr3 molecules into the trigonal prismatic 
voids at (1/3,2/3,1/2) and (2/3,1/3,1/2) [Fig. [f(b)]. Based on 
the modeled crystal structure they calculated an X-ray diffrac- 
tion pattern almost identical to the observed one. The crystal 
structures can be considered as a sequence of alternating lay- 
ers of Cgo and haloform molecules. As an example, Fig. [jjd) 
shows the layer structure of Cgo ■ 2 CHCI3. 

Very recently, Dinnebier et al. ^ performed high resolu- 
tion X-ray powder diffraction analyses at various tempera- 
tures and determined the crystal lattices of different phases 
of C60 • 2 CHX3 (X=Cl,Br). In agreement with Jansen and 
Waidmann and Collins et al., the room temperature crystal 
lattice obtained is hexagonal and has P6/mmm space group 
symmetry. At T w 200 K the haloform intercalated fullerenes 
undergo a first-order phase transition towards a monoclinic 
phase of C2/m symmetry which has the double volume of the 
room temperature phase. The alternating layers of haloform 
and C60 molecules are shifted with respect to each other such 
that the molecules are displaced from their equilibrium posi- 
tions in the room temperature phase. Below ~150 K a second 
first-order phase transition occurs and the crystal structures of 
Cgo ■ 2 CHX3 are triclinic and have PI symmetry [[|[. This 
low temperature (LT) phase has similar cell dimensions as the 
hexagonal room temperature structure (Table Q). Throughout 
this paper this structure is labeled with LT. 

A comparison of the hexagonal unit cells of C60 • 2 CHX3 
(X=Cl,Br) [Fig. |l](a) and |](b)] with the primitive unit cell 
of pure fee Cgo [Fig. [l|(c)] emphasizes the structural changes 
caused by the intercalation with haloform molecules. The cu- 
bic closed-packed (ccp) lattice of Cgo molecules transforms to 
a hexagonal close-packed lattice of C60 molecules. As a result 
of the intercalation the length of the primitive cell is increased 
by 1 to 4 %. 



III. COMPUTATIONAL METHOD 

All calculations are performed with the DMol 3 approach 
Jig 1 , jig ]. The calculations make use of the gradient cor- 
rected density functional of Becke [|l7|] for the exchange con- 
tribution, combined with the Perdew Wang 1991 functional 
JT^ ] for the correlation contribution. All electrons are de- 
scribed explicitely with a double numeric basis set augmented 
with polarization functions (DNP): 2slp for H, 3s2pld for C, 
4s4p2d for CI and 5s4s2d for Br. The cutoff radius to gen- 
erate the atomic basis functions is set to 7.0 bohr (1 bohr = 
5.29177 x 10~ n m). The electron densities and the energies 
are calculated self consistently, latter within an accuracy of 
10 -5 hartree (1 hartree= 2625.5 kJ mol -1 ). The numerical in- 
tegration is performed on medium sized meshes. For the Bril- 
louin zone integration irreducible k vector meshes of different 
sizes are used Jig 1 ]. Reliable energies and atomic positions 
are obtained with a k vector mesh of (n x ,n y ,n z ) = (4,4,2), 
where n, denotes the increments along each of the primitive 
reciprocal lattice vectors, G; (Fig. ||) [jlq]. This k vector mesh 
amounts to 20 symmetry-inequivalent sampling points. The 
calculations of the DOS required a finer k vector mesh of 
(n x ,n y ,n z ) = (8,8,4), i.e. 132 irreducible k points. Calcula- 
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FIG. 1: Calculated configurations of (a) CHCI3 and (b) CHBr3 molecules intercalated in hexagonal lattices of Cm molecules. Both crystal 
structures have PI symmetry and are based on the X-ray powder diffraction analysis of Jansen and Waidmann mn. (c) The molecules in 
pure Cso form an fee lattice and the crystal structure has Fm3 symmetry. Depicted is the primitive unit cell of fee Cso- (d) The structures of 
hex. Cgo • 2 CHCI3 and hex. C60 ■ 2 CHBi'3 can be considered as a sequence of alternating layers of C60 molecules and haloform molecules 
perpendicular to the c vector. 



tions on fee C60 on the other hand, require only 35 symmetry- 
unique sampling points, (n x ,n y ,n z ) = (8,8,8). Using the self- 
consistent electronic densities the electronic band structures 
of fee C 6 o and C 6 o • 2 CHX 3 (X=Cl,Br) are calculated at 59 
and 153 reciprocal lattice vectors, respectively, at the surface 
of the irreducible wedge of the appropriate Brillouin zones. 

Doping with holes (electrons) is simulated by simply re- 



moving (adding) an appropriate number electrons from the 
unit cell. For instance, to simulate a doping of three holes per 
C60 molecule, three electrons are removed and the total charge 
of the unit cell is denoted Q = +3. The resulting positive or 
negative charge of the system is compensated by a background 
counter charge (jellium approach). Considered are doping lev- 
els of Q = 0, Q = ±1 and Q = ±3. The charging procedure 



4 



forms systems with open electron shells. Therefore, all calcu- 
lations on the charged systems are performed within the un- 
restricted Kohn-Sham formalism to account for possible spin 
polarization effects. 

The crystal structures introduced in the next section for 
which the electronic properties are determined, are optimized 
by keeping the cell parameters fixed and allowing the atoms 
to move. This concerns neutral as well as charged sys- 
tems. The atomic positions are optimized employing delo- 
calized internal coordinates Jl9[]. The convergence criteria of 
the optimization procedures are 10~ 5 hartree for the energy, 
3 x 10~ 4 hartree/bohr for the maximum energy gradient and 
3 x 10~ 3 bohr for the maximum atomic displacement. 



IV. RESULTS AND DISCUSSION 
A. Optimized crystal structures 

Calculations are performed on the hexagonal HT structures 
of Jansen and Waidmann [[!]] as well as on the triclinic LT 
structures of Dinnebier et al. Because the atomic posi- 
tions in the unit cells of the HT structures are unknown, sev- 
eral different atomic configurations are designed as starting 
points for the structure optimizations. Physically and chem- 
ically reasonable configurations can only be designed when 
the experimental determined space group symmetry P6/mmm 
is reduced. 

One designed configuration of Ceo ■ 2 CHCb has P3 sym- 
metry. The two chloroform molecules in the unit cell are 
located in the trigonal prismatic voids at (1/3,2/3,1/2) and 
(2/3,1/3,1/2), respectively, and are related by inversion sym- 
metry. One six-fold axis of the Cgo and the C3-axes of the 
CHCI3 molecules coincide with the c axis of the hexagonal 
lattice. The two C-H bonds point towards the — c and +c lat- 
tice vector, respectively. Three other start configurations are 
derived from the P3ht configuration by exchanging the hy- 
drogen atoms of the CHCI3 molecules with arbritrary chlo- 
rine atoms. Two of these configurations have PI symmetry 
and one has PI symmetry. Due to the exchange of H and CI 
atoms the CHCI3 molecules are differently oriented and prob- 
ably induce also an orientational change of the C60 molecules. 
Throughout this article all configurations are named by their 
symmetry labels, i.e. the configurations formed are termed 
PIht, Plffr> P3ffT an d Pl#r- The most stable configuration 
of C60 • 2 CHCI3, i.e. the one with the lowest calculated en- 
ergy, is used to form the starting configuration for the structure 
optimizations of Cgo ■ 2 CHBr3 . 

For the calculations on pure fee C60 the crystal structure of 
David is employed [jljj]. The space group of this structure is 
Fm3, i.e. all molecules have the same orientation. Its primitive 
unit cell contains one molecule [Fig. |l](c)] and has a length of 
a = 9.949 A (Table §). 

It is worth noting that both CHCI3 and CHBr3 have per- 
manent electric dipole moments of /j ~ 1 D [gG] and the C60 
molecules are highly polarizable (a ~ 80 A 3 [pip). The charge 
distribution within the haloform molecules is such that the 
halogen atoms are charged partially negatively and the hydro- 



gen atoms are charged partially positively. Hence, the sta- 
bility of C60 • 2 CHX3 (X=Cl,Br) is presumably mainly due 
to dipole-induced dipole interactions rather than by van-der- 
Waals interactions between the molecules. 

Table [n| summarizes the relative energies, AE, between the 
optimized crystal structures of neutral C60 ■ 2 CHCI3 and 
C60 ■ 2 CHBr3. All energies are with respect to the most 
stable configurations which have relative energies of AE= 
kJ-mol -1 . The most stable structures of C60 ■ 2 CHCI3 are 
the PIht configurations [Fig. |](a)]. The energy difference 
to the next stable configuration (P3ht) is only 3 kJ-mol -1 . 
Although the LT configurations correspond to temperatures 
below ^150 K and the HT structures to room temperature, 
the calculations on C60 ■ 2 CHCI3 yields a lower energy 
for the PIht configuration than for the Pl/r configuration. 
The energy difference is A£= 28 kJ-mol -1 . In the case of 
C60 ■ 2 CHBr3 both structures, PIht an d Pl.tr> have nearly 
identical energies. 

These energy considerations are made for zero temperature 
and, therefore, do not include any entropic effects. The en- 
ergies shall give only a clue about the shape of the potential 
energy surface of C60 ■ 2 CHX3 (X=Cl,Br). The small energy 
differences between the configurations give rise to a shallow 
potential energy surface provided that the energy barriers of 
molecular re-orientations are of the same order as the rela- 
tive energies. Solid-state NMR measurements predicted for 
the re-orientations of the C60 molecules in C60 • 2 CHBr3 ac- 
tivation energies of 6 kJ-mol _1 |Q]. Nevertheless, in order to 
give a more precise answer of the relative stability between 
the configurations Cgo ■ 2 CHCI3 and C60 • 2 CHBr3 at least 
the zero-point vibrational energy and the vibrational partition 
function need to be calculated. 

All optimized structures of C60 ■ 2 CHX 3 (X=Cl,Br), 
obtained from the HT and LT data have almost identical 
intramolecular atomic distances (in pm): R(C-C) = 145, 
R(C=C) = 140, R(C-U) = 109, fl(C-Cl) = 177-179, K(C-Br) 
= 194-196. These distances are equivalent to the appropriate 
distances of free molecules in the gas phase. 

The main differences between the calculated structures can 
be found when comparing the intermolecular distances. The 
differently oriented haloform molecules in the start config- 
uration of the optimization procedure induce changes in the 
orientation of the fullerene molecules. We note that for both 
chloroform and bromoform intercalation the high resolution 
X-ray measurements [|^] and the calculated LT structures (this 
work) are almost identical. The largest differences between 
these structures are in the intermolecular distances and are 
smaller than 5 pm. This means, the optimization of the atomic 
positions of the LT structures with fixed lattice parameters 
does not cause any drastical structural changes. 

Table [n] summarizes selected intermolecular distances of 
the optimized crystal structures of neutral fee Cgo and of neu- 
tral C 60 • 2 CHX 3 (X=Cl,Br). While the expansion of the lat- 
tice of C60 molecules caused by the intercalation is accompa- 
nied by an increase of the unit cell length of only 1 to 4 % the 
changes of the shortest distances between Cgo molecules vary 
between 1 and 18 %. This indicates that there are changes 
in the orientation of the molecules. Due to the larger effec- 
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TABLE II: Relative energies between the different configurations (A£, in kj-mol ) and selected intermolecular atomic distances (in pm) 
of several calculated structures of neutral fee Cgo and C50 ■ 2 CHX3(X=Cl,Br). Summarized are the shortest distances and, if given, the 
corresponding average values in parenthesis. Note, the distances of fee C(,q are related to the primitive unit cell [Fig. Il|(c)] 





AE 


C- • C (Cm) 
along b 


along a-\- b 


along c 


C--X a > 


X- -X b ) 


C -H b) 


fee CtA c) 
ice 
















Fm3 




299 301 f ' 


315 f ) 


342 f ' 








C 60 ■ 2 CHCI3 
















Pl H r d) 





340, 302 


340 


354 


301 (339) 


308 (348) 


262 (347) 


P3ht d) 


3 


306, 306 


306 


362 


355 


350 


334 (359) 


Pi' d ) 


11 


304, 302 


346 


355 


317 (338) 


298 (328) 


295 (345) 


Pl HT d ' 


14 


303, 302 


351 


356 


319 (342) 


305 (326) 


303 (357) 


Pl Lr e ) 


28 


326, 324 8) 


340 8) 


3618) 


339 (347) 




269 (335) 


C 60 ■ 2 CHBr 3 
















Pl H r d) 





354, 335 


343 


360 


300 (341) 


318 (342) 


251 (342) 


Pl Lr e > 


1 


332, 334 s) 


354 8) 


385 8) 


343 (349) 




286 (337) 



') Distances shorter than 357 pm. b ) Distances shorter than 400 pm. c ' Crystal structure of David p4[]. d ) Crystal structure of Jansen and 
Waidmann [jlj]. e ^ Crystal structure of Dinnebier et al. [g]. f ' The c vector of the hexagonal unit cell corresponds to the b vector of the 
primitive fee unit cell. g ) The c vector of the hexagonal unit cell corresponds to the b vector of the triclinic unit cell used in Ref. pp. 



tive volume of CHB1-3 the shortest interatomic distances be- 
tween neighboring Cgo molecules are on the average larger 
for C 6 o • 2 CHBr 3 than for C 6 o ■ 2 CHCI3 and fee C 6 o- But 
the largest changes are found for the shortest distance between 
C atoms of nearest neighboring C60 molecules along the cell 
vector a, R(C---C) a , of the PIht structures, which are 41 
(X=C1) and 53 pm (X=Br). Note, that the distances of fee 
Cgo are related to the primitive unit cell [Fig. [j](c)]. 

All the shortest C- ■ ■ C distances along the c vector, 
R(C---C) C , of the HT as well as the LT structures are 
larger than the shortest C- • C distances in the (0001) plane, 
typical for layered system structures. By far the largest 
R(C- ■ ■C) e value of 385 pm is for the Plw configuration of 
C 6 o ■ 2 CHBr 3 . 

Both in fee Cgo and the HT structures of C60 • 2 CHCI3 
there is at least one value of R(C- ■ - C) a or R(C- ■ -C)b which 
is around 300 pm. These distances are about 40 pm smaller 
than the sum of the van-der-Waals radii of two carbon atoms 
(Rvdw(C)= 170 pm J22]]). This is evidence for an overlap 
between 71 orbitals of neighboring Ceo molecules. While 
the Pl/T and most of the HT structures of Ceo ■ 2 CHCI3 
have almost identical values for R(C---C) a and R(C---C)t, 
only PIht Ceo • 2 CHCI3 has two different values, where 
R(C- ■ -C) a is about 40 pm larger than R(C- ■ -C)&. Because 
of symmetry restrictions for P3 ht C(,q ■ 2 CHCI3, all the 
shortest intermolecular C - ■ C distances in the (0001) plane, 
R(C- ■ ■ C%, R(C- - C) h and R(C- ■ ■ C) a+b , are equivalent. 

On the other hand the LT configuration of C60 • 2 CHCI3 
(PlLr) does not have any interatomic distance between neigh- 
boring C60 molecules which is shorter than 324 pm. But the 
values of R(C- ■ ■ C)„ and R(C- ■ ■ C)b are also almost identical, 
similar to fee Ceo- 

In both structures of Cgo • 2 CHBr 3 considered only 
R(C- ■ -C)b are shorter than 340 pm, the sum of the van-der- 
Waals radii of two carbon atoms. Similar to Cgo ■ 2 CHCI3 the 
PIht structure of C60 • 2 CHBr3 has two different values for 



R(C- ■ -C) a and R(C- ■ -C)b while for PIlt these distances are 
almost identical. 

For all configurations both the average values and the small- 
est values of the distances between the halogen atoms X 
(X=Cl,Br) and the atoms of the Ceo molecules are smaller than 
the sum of the corresponding van-der-Waals radii {R M iw{CY)= 
175 pm, R vc iw(Bi)= 185 pm [Q). The PIht configurations 
of C60 • 2 CHCI3 and Cgo • 2 CHBr3 have the smallest value of 
R(C- ■ -X)f« 300 pm. This indicates a significant intermolecu- 
lar overlap between orbitals of Cgo and CHX3. 

Due to their partial negative charge the halogen atoms 
of neighboring haloform molecules should avoid each other. 
Therefore, it is peculiar that the majority of the configurations 
have halogen atoms X that are closer than twice of R v <jw(X). 
Only the PIlt structures do not have any X---X distance 
smaller than 400 pm. 

All configurations have similar average values of the non- 
bonding C- ■ ■ H distances. The PIht structures have the short- 
est non-bonding C - H distances. Apparently, the partially 
positively charged H atoms are attracted by the negatively 
charged 7C electron system of the Cgo molecules. 

The electronic structure calculations of the following sec- 
tions are performed representatively for the most stable 
(Pl//r) and most unstable (PIlt) configurations as deter- 
mined in this section. It is assumed that the other structures 
have similar properties which can be derived or approximated 
from the properties of the PIht and PIlt configurations. 

B. Electronic structures 

1. Band structures and DOS 

Fig. [| shows the valence and conduction bands of neutral 
PIht C 60 • 2 CHX 3 (X=Cl,Br). The corresponding Bril- 
louin zone is depicted in Fig. The systems have indirect 
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FIG. 2: Hexagonal Brillouin zone with reduced symmetry. 



band gaps of 1.2 (X=C1) and 1.3 eV (X=Br). Since neutral 
Ceo • 2 CHX3 has empty conduction bands the systems are 
non-metallic. Their valence bands (VB) and conduction bands 
(CB) consist of five and three narrow bands, respectively, in 
accordance with the respective degeneracy of the constitutive 
molecular orbitals. The total DOS of the CB clearly show 
three nearly separated peaks, each belonging to one of the 
three bands. None of the VB and CB are degenerate at any 
k vector of the band structure. The largest band dispersions 
are along k x and k y . Table |Th| collects the width of the VB and 
CB and the smallest indirect gaps between VB and CB. The 
widths of the VB and CB of PI HT C 60 • 2 CHBr 3 are 0.4 and 
0.3 eV, respectively. The VB of Pl„ r C 60 ■ 2 CHCI3 is 0.1 eV 
broader, because the fullerene molecules are closer to each 
other (see Tab. ||) and, hence, the overlap of the tc orbitals of 
neighboring molecules is larger. Note, that the band structure 
of neutral Cgo ■ 2 CHBr 3 has a distinguishing feature namely 



that a single narrow band appears just 0.4 eV below the VB. 
In the case of C60 • 2 CHCI3 the electronic states below the 
VB are rather a group of bands. The gap between the latter 
and the VB is about 0.5 eV. 

The LT structures have similar band structures as the PIht 
structures but with band gaps of 1.1 (C60 ■ 2 CHCI3) and 1.2 
(C 60 ■ 2 CHBr 3 ). For both, C 60 ■ 2 CHCI3 and C 60 ■ 2 CHBr 3 , 
the width of the VB is 0. 1 eV larger than that of the appro- 
priate HT structure in agreement with the intermolecular dis- 
tances of Table Q. The widths of the CB and the gaps between 
the VB and the upper edge of the next lower lying electronic 
states are identical for the HT and LT structures. 

For fee C60 the calculated band gap is only 1 .0 eV and the 
valence and conduction bands are 0.1 to 0.3 eV broader than 
the VB and CB of the HT and LT structures. Both, the smaller 
band gap and the broader bands are due to smaller distances 
between the fullerene molecules as compared to C60 ■ 2 CHX3 
(X=CLBr). 

In order to disentangle the contribution of the chloroform 
molecules to the electronic structure of the C60 subsystem in 
Cgo • 2 CHCI3, calculations are also performed on a hypo- 
thetical hexagonal lattice of Cgo molecules in the geometry of 
PIht Ceo ■ 2 CHCI3. This model system is designed by sim- 
ply removing the two CHCI3 molecules from the unit cell of 
PIht Ceo • 2 CHCI3. The shape and the positions of the va- 
lence and conduction bands almost do not change. The band 
gap of 1 .2 e V and the widths of the VB and CB are identical 
to those of PIht C(,q ■ 2 CHCI3. Hence, the chloroform as 
well as the bromoform molecules apparently do not affect the 
band structure of the materials substantially. 

One of the well-known drawbacks of DF calculations is the 
underestimation of band gaps. For instance, the experimental 
determined band gap of fee Cgo is 2.3 eV [ ^3| ] whereas the 
calculated gap is less than half of that (vide supra). Hence, 
the band gaps presented in this article are only lower bounds. 
The problem of too small band gaps in DF calculations can be 
cured according to Slater transition-state methods g5j]. 



TABLE III: Comparison of width and smallest indirect gaps for the 
valence and conduction bands, VB and CB, respectively (in eV) cal- 
culated for fee C 6 o and neutral C 6 o • 2 CHX 3 (X=Cl,Br). 





band gap 
VB - CB 


band width 
VB 


CB 


fee C 60 a ) 


1.0 


0.7 


0.6 


C 60 ■ 2 CHCI3 








PIht b) 


1.2 


0.5 


0.4 


Pl Hr , without CHCI3 b )< c ) 


1.2 


0.5 


0.4 


PIlt- d) 


1.1 


0.6 


0.4 


C 60 • 2 CHBr 3 








PlflT 


1.3 


0.4 


0.3 


PIlt- d) 


1.2 


0.5 


0.3 



a ) Crystal structure of David jl4j|. b ' Crystal structure of Jansen and 
Waidmann [Q]. c ' The two CHCI3 molecules are simply deleted 
from the unit cell of the PIht structure and the remaining Cgo 
molecule is re-optimized. d ' Crystal structure of Dinnebier etal. 



2. Partial DOS 

Fig. |] shows the partial density of states (PDOS) for spe- 
cific atoms of PIht C 60 ■ 2 CHX 3 (X=Cl,Br). The PDOS 
is the decomposition of the electron bands into contributions 
of atomic orbitals (AO). The atoms of C60 and the halo- 
form molecules chosen for the PDOS have the smallest in- 
termolecular C- • X distances. As expected, the VB and the 
CB of PIht C 6 o ■ 2 CHCI3 and C 6 o • 2 CHB13 mainly con- 
sist of carbon 2p n AO. Surprisingly, also chlorine 3p and 
bromine 4-p AO contribute significantly to the PDOS of PIht 
C60 ■ 2 CHCI3 and C60 • 2 CHBr3, respectively. However, 
the latter two contributions are smaller than that of the C60 
atoms. The single electronic band just below the VB of PIht 
Ceo ■ 2 CHBr3 (Fig. ||) mainly consists of bromine Ap AO. 

The calculated PDOS for the LT structures leads to con- 
clusions similar to those for the HT structures. However, 
the contributions of the p-type orbitals of the halogen atoms 
are smaller and the contribution of fullerenes correspondingly 



--60 



2 CHCh 



-60 



2 CHBr. 




wavevector DOS wavevector DOS 

FIG. 3: One-electron band structures and total density of states (DOS) of neutral PIht Cgo ■ 2 CHX3 (X=Cl,Br). All energies are with respect 
to the Fermi energy, Ef . The band structures depict the three conduction bands, the five valence bands and the next state below the valence 
bands. For clarity, the bands are gray and black in alternating order. 



larger, because the LT structures have larger intermolecular 
C- • X (X=Cl,Br) distances on average. 

Consequently, the haloform molecules do not only have the 
function of inert spacers to expand the Cgo lattice as inferred 
in Section 



IV Bl 



but might also support charge transfer and 



electron-phonon coupling. 



3. Charge distribution 

This section concerns only neutral PIht Cgo ■ 2 CHCI3 
because neutral PIht C60 • 2 CHBr3 and neutral PIlt 
Cgo • 2 CHX3 (X=Cl,Br) have akin charge density distribu- 
tions. 

Figs. |5](a) and^(b) show contourplots of the valence charge 
density in the (1100) and (1 100) planes, respectively. The va- 
lence charge density is the sum of the squares of the five high- 
est occupied crystal orbitals of neutral Cgo ■ 2 CHCI3. The 
integral over the Brillouin zone is obtained using the mid- 
point approximation (r point). In the corner of the maps are 
parts of the C60 molecules represented by their Jt electron sys- 
tems. Clearly visible are the carbon 2p K AO. In the center of 
Fig^(a) valence charge density that belongs to the two CHCI3 
molecules can be seen. Both density plots show considerable 
overlap of the p K orbitals of adjacent C60 molecules as well as 
overlap of the chlorine p orbitals with the 71 electron systems 
of the C60 molecules. The latter overlap is smaller than the 
overlap between the Jt electron systems, yet relevant. 

The contour plots of Fig. || show the molecular deforma- 
tion density, Pdiff, in particular planes. The molecular defor- 
mation density is defined as the difference between the total 
electron density of C60 ■ 2 CHCI3 and the sum of the densities 



of the subsystems C60 and CHCI3. Note, that both subsystems 
have the same geometry as in C60 • 2 CHCI3. The molecu- 
lar deformation density discloses changes of the electron den- 
sity of the C60 and chloroform molecules due to intercalation, 
e.g. polarization effects or charge transfer from one molecule 
to another. 

Fig. ^(a) shows the changes of the electron density induced 
by chlorine atoms that have smallest interatomic distances to 
the C60 molecules. One chlorine atom pulls electron density 
away from C60 molecule and the other chlorine atom pushes 
electron density towards a Ceo molecule. Both, Fig. ^(b) and 
Fig. ||(c), depict the polarization of carbon-carbon Jt bonds 
due to the two partially positively charged hydrogen atoms of 
the chloroform molecules. These hydrogen atoms have also 
the smallest interatomic distances to Cgo atoms. The electron 
density surrounding the 7C-bonds is increased at the expense of 
the density around the hydrogen atoms. 

Note that the deformation density is of the order of 10~ 4 
atomic units and, hence, very small. However, the contour 
plots of Fig. U clearly validate the presumption of dipole- 
induced dipole interactions between the haloform and the Cgo 
molecules. 



C. Changes of the electronic structures due to charging 

In order to clarify how electron and hole doping affect the 
electronic structure and perhaps also the atomic configura- 
tions of haloform intercalated Cgo, calculations are performed 
on the charged PIht and PIlt structures of Cgo ■ 2 CHX3 
(X=Cl,Br) as well as on charged fee Cgo- The structural 
changes of the PIht and PIlt configurations of Cgo ■ 2 CHX3 
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FIG. 4: Partial DOS of the bands depicted in Fig. Q. All energies 
are with respect to the Fermi energy, Ep . The partial DOS show the 
contributions of the s-, p- and d-type of orbitals of specific atoms of 
the C60 molecule, C, and of the haloform molecules, CI and Br, to 
the total DOS of P1 HT C 60 ■ 2 CHX 3 (X=Cl,Br). Note the different 
ranges of the ordinates. 



due to charging are found to be small. The charging has no ef- 
fect on the intramolecular distances. For example Jahn-Teller 
distortion could not be found. The changes in intermolecular 
distances of Table [y| are of the order of 2 pm and are caused 
by small changes in the rotational orientation of the haloform 
and C60 molecules. 

The Hirshfeld population analysis p6| ] predicted that in 




FIG. 5: Contour plots of the sum of the squares of the five highest 
occupied crystal orbitals at the T point of neutral PI//7- Cgo • 2 CHCI3 
in the [(a)] (1 100) and [(b)] ( 1 100) planes. The iso-lines depict values 
from 10~ 5 to 5.12 x 10~ 3 each increased by a factor of two (1 a.u. = 
6.7483 x 10 30 electrons ■ m~ 3 ). The darker the gray the higher is the 
electron density. 



each case 80 % to 90 % of the charge is located on the Ca) 
molecules and the remaining part on the haloform molecules. 

Fig. ^ shows the total DOS of the five valence and three 
conduction bands of neutral Cgo • 2 CHX3 (X=Cl,Br) and their 
modifications upon hole and electron doping. In the follow- 
ing the terms VB and the CB are related to the corresponding 
states of the neutral systems. Overall the DOS of the VB and 
CB are narrow peaks with steep flanks. However, for the same 
charge the shape of the DOS of the HT and LT structures show 
differences. For almost all charges the CB has three distinct 
peaks which indicates that it consists of three narrow bands. 
The widths of the CB are in the range of 0.3-0.5 eV with- 
out any clear trend with respect to charging. The width of the 
DOS of the five valence bands is between 0.4 to 0.6 eV and 
the single bands are not resolved. 

Significant are the shifts in energy of both VB and CB. 
Upon electron doping both the VB and the CB are shifted to 
higher energies whereas hole doping causes a shift to lower 
energies. The gap between the VB and the CB has its maxi- 
mum at zero charge and decreases for positive as well as for 
negative charging by 0.1-0.2 eV. The gap between the VB 
and the states lower in energy as well as the gap between 




FIG. 6: Contour plots of the molecular deformation density, pdiff' plotted for relevant planes of PIht C^q ■ 2 CHCI3 (see text for the 
definition of Pdiff)- ( a ) shows the changes of the total density induced by CI atoms, and (b) and (c) induced by the two H atoms in the unit cell. 
The atoms in and near the plane are marked by open circles. The iso-lines represent values from — 10~ 4 to 10 -4 a.u. in steps of 10~ 5 a.u. . 
Dashed lines are for negative and solid lines for positive values. The light-gray solid lines sketch parts of the CHCI3 and C50 molecules. 
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FIG. 7: Total density of states (DOS) around the Fermi energy, Ef, of the neutral and charged PI//7 and Pl/j structures of C50 ■ 2 CHX3 
(X=Cl,Br). 2 is the total charge per unit cell. The vertical solid lines represent Ef. The corresponding total DOS at Ef are summarized in 



Tab.Q 



the CB and the states higher in energy are also changed. In 
most of the cases the gap between the CB and the next higher 
bands is larger than 0.6 eV. However, for both structures of 
C60 • 2 CHBr3 and Q = — 3 this gap is only 0.2 eV. For each 
structure of C60 ■ 2 CHX 3 (X=Cl,Br) the gap between the VB 
and the next lower states is almost constant upon electron dop- 
ing whereas for hole doping the VB are shifted towards lower 



energy and the next lower states towards higher energy. For 
Plffr C60 • 2 CHBr 3 and Q = +3 the former single band below 
the VB has merged with the latter, i.e. the gap has vanished. 
The VB is formed of six instead of five states. However, there 
is still a small gap of ^0.1 eV between these six bands and 
the next lower states. More drastic are the changes in the case 
of PI LT C 60 • 2 CHBr 3 and Q = +3 where the VB and the 



10 



TABLE IV: Total density of states (DOS) at the Fermi energy of 
charged pure fee C^q and C^q ■ 2 CHX3 (X=Cl,Br). All values are in 
states per unit cell and eV. 
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a ' Total charge per unit cell. b ) Crystal structure of David [Q], 
c ) Crystal structure of Jansen and Waidmann [|l|]. d ' The two CHCI3 

molecules are simply deleted from the unit cell of the PIht 
structure and the remaining Cgnmolecule is re-optimized. e ' Crystal 

structure of Dinnebier et al. ph. f ' Average of PIht and Plrj-. 



next lower states form a continuum. 



Table IV summarizes the calculated total DOS at the Fermi 
energy of charged pure and haloform intercalated C60 for Q = 
±1 and Q = ±3. The accuracy of the calculated DOS is about 
± 1 state per unit cell and e V. 

There is a clear trend that negatively charged systems have 
smaller DOS than positively charged systems. The small dif- 
ferences in shapes of the DOS of the PIht and PI u structures 
(Fig. become apparent in remarkable differences between 
DOS at the Fermi energy. The DOS at Ep of the PI ht struc- 
tures are for some charges larger and for other charges smaller 
than that of the PI £7- structures. The discrepancies are up to 
12 states per eV and unit cell and are primarily a result of 
different intermolecular distances (Table caused by differ- 
ently oriented haloform and Cgo molecules. Furthermore, due 
to narrow VB and CB peaks the DOS is sensitive with respect 
to the position of the Fermi energy. 

However, at finite temperature the PIht as well as the Pl/,7 
configuration can occur and the DOS at the Fermi energy for a 
particular charging would be an average of both Except 
for Q = +1 the averaged DOS is increasing as DOS(fcc Cm) 
< DOS(C 60 • 2 CHCI3) < DOS(C 60 • 2 CHBr 3 ) (Table g). 
But also for Q = +1 Ceo ■ 2 CHCI3 has on average a smaller 
DOS than C 60 -2CHBr 3 . 

Table ^ includes also the DOS at the Fermi energy of the 
hypothetical hexagonal lattice of pure 0*6 in the geometry of 
PIht C 60 ■ 2 CHC1 3 (see Section fVBl\). For all charges the 



calculations yields DOS that are very similar to that of PIht 
Cgo ■ 2 CHCI3. Hence, the orbitals of the CHCI3 molecules 
which interact with the K system of the Cgo subsystem have a 
negligibly influence on the DOS of the entire charged system. 



V. SUMMARY AND CONCLUSIONS 

Using density functional methods we have determined the 
electronic properties (band structure, partial and total density 
of states, charge distribution) of pure C60 and C60 • 2 CHX3 
(X=Cl,Br) after optimization of their atomic structures. The 
changes resulting from doping with one and three elec- 
trons/holes per unit cell have also been examined. 

The calculations yielded for un-doped chloroform and bro- 
moform intercalated C(,o, C(,o • 2 CHX3 (X=Cl,Br) indirect 
gaps between the valence bands (VB) and conduction bands 
(CB) larger than 1 eV. Both C 60 ■ 2 CHC1 3 and C 60 • 2 CHBr 3 
are narrow band materials, i.e. most of the electronic states are 
mainly localized on the molecules. The calculated width of 
the VB and CB are 0.4-0.6 eV and 0.3-0.4 eV, respectively. 
The orbitals of the haloform molecules have a considerable 
overlap between the Jt orbitals of the fullerene molecules and 
the p-type orbitals of halogen atoms significantly contribute 
to the VB and CB of Cgo • 2 CHX3. Both compounds are 
mainly stabilized by dipole-induced dipole interactions rather 
than by van-der-Waals interactions. The intrinsically differ- 
ently charged atoms of the haloform molecules cause a polar- 
ization of the fullerene molecules. 

Doping with charge carriers turns the C60 • 2 CHX3 
(X=Cl,Br) to metals. As a result of the intermolecular overlap 
of the CHX 3 and C 60 orbitals, 10 to 20 % of the charge of 
the doped systems is on the haloform molecules instead of be- 
ing completely localized on the fullerene molecules. Charging 
with electrons leads to a shift of both VB and CB to higher 
energies whereas the doping with holes shift them to lower 
energies. The corresponding band widths undergo relatively 
small changes. At a charging of Q = +3 per C60 • 2 CHX3 the 
gap between the VB and the states just below the VB is signif- 
icantly reduced. In the case of Cgo ■ 2 CHBr3 the VB has even 
merged with states lower in energy. The calculated density of 
states at the Fermi energy for different chargings are clearly 
smaller for electron than for hole doping and are larger for 
C 6u • 2 CHBr 3 than for C 60 • 2 CHCI3. Calculations on dif- 
ferent crystal structures of C(,o • 2 CHCI3 and C(,o • 2 CHBr3 
revealed that the density of states at the Fermi energy are sen- 
sitive to the orientation of the haloform and C60 molecules. 
At a given charge the differences between the various crystal 
structures are up to 12 states per eV and unit cell. 

At a charging of Q = +3, which resembles the supercon- 
ducting phase of pure Cgo and C(,o • 2 CHX3 (X=Cl,Br), 
we calculated DOS at the Fermi energy that increases as 
DOS(Ceo) < DOS(C 60 • 2 CHCI3) < DOS(C 60 • 2 CHBr 3 ). 
On the other hand, the tight-binding calculations of Dinnebier 
et al. |Q| yielded DOS that have the reverse order. This dis- 
crepancy is related to the different methods considered. The 
method used in the present work is more accurate than the 
tight-binding approach. Nevertheless, our calculations clearly 
supports the conclusion of Dinnebier et al. (§] that the DOS 
alone cannot account for the observed increase of T c upon 
CHCI3 intercalation and substitution of CI by Br. According 
to a different scenario proposed by Bill and Kresin [ |28| ] the 
increase of T c upon intercalation is caused by an additional 
contribution to the superconducting pairing interaction which 
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comes from the coupling of charge carriers with the vibra- 
tional manifold of the intercalated molecules. The difference 
in vibrational spectra of CHCI3 and CHBr3 molecules [ p9[ ] 
leads to a noticeable increase of T c upon Cl^Br substitution. 

The present calculations have all been performed on a 
three-dimensional crystal in absence of an applied electric 
field. However, the presence of an electric field may lead to a 
re-orientation and polarization of the haloform and fullerene 
molecules. Such structural changes affect the electronic struc- 
ture of the system and may influence vibronic couplings as 
well. It would be of interest to complement the study pre- 



sented in this article with calculations on C60 ■ 2 CHX3 
(X=Cl,Br) in an electric field. A similar problem has been 
examined for pure C^o [pOj] . 
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